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Abstract — A  hypergraph  is  a  set  V  of  vertices  and  a  set 
of  non-empty  subsets  of  V,  called  hyperedges.  Unlike  graphs, 
hypergraphs  can  capture  higher-order  interactions  in  social  and 
communication  networks  that  go  beyond  a  simple  union  of 
pairwise  relationships.  In  this  paper,  we  consider  the  shortest 
path  problem  in  hypergraphs.  We  develop  two  algorithms  for 
finding  and  maintaining  the  shortest  hyperpaths  in  a  dynamic 
network  with  both  weight  and  topological  changes.  These  two 
algorithms  are  the  first  addressing  the  fully  dynamic  shortest 
path  problem  in  a  general  hypergraph.  They  complement  each 
other  by  partitioning  the  application  space  based  on  the  nature  of 
the  change  dynamics  and  the  type  of  the  hypergraph.  We  analyze 
the  time  complexity  of  the  proposed  algorithms  and  perform 
simulation  experiments  for  both  random  geometric  hypergraphs 
and  the  Enron  email  data  set.  The  latter  illustrates  the  application 
of  the  proposed  algorithms  in  social  networks  for  identifying  the 
most  important  actor  based  on  the  closeness  centrality  metric. 

I.  Introduction 

A  graph  is  a  basic  mathematical  abstraction  for  modeling 
networks,  in  which  nodes  are  represented  by  vertices  and  pair¬ 
wise  relationships  are  represented  by  edges  between  vertices. 
A  graph  is  thus  given  by  a  vertex  set  V  and  an  edge  set  E  con¬ 
sisting  of  cardinality-2  subsets  of  V.  A  hypergraph  is  a  natural 
extension  of  a  graph  obtained  by  removing  the  constraint  on 
the  cardinality  of  an  edge:  any  non-empty  subset  of  V  can  be 
an  element  (a  hyperedge)  of  the  edge  set  E  (see  Fig  1).  It  thus 
captures  group  behaviors  and  higher-dimensional  relationships 
in  complex  networks  that  are  more  than  a  simple  union  of 
pairwise  relationships.  Examples  include  communities  and 
collaboration  teams  in  social  networks,  document  clusters 
in  information  networks,  and  cliques,  neighborhoods,  and 
multicast  groups  in  communication  networks. 


Fig.  1.  An  example  hypergraph  with  4  hyperedges:  (vi,  V2,  V3,  vg), 
(V2,V3,V4,V5),  (t>6,V7,t;8,l’9),  and  (v5,V8,Vg). 

While  the  concept  of  hypergraph  has  been  around  since 
1920’s  (see,  for  example,  [1]),  many  well-solved  algorithmic 
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problems  in  graph  theory  remain  largely  open  under  this  more 
general  model.  In  this  paper,  we  address  the  shortest  path 
problem  in  hypergraphs. 

A.  Shortest  Path  Problem  in  Graphs 

The  shortest  path  problem  is  perhaps  one  of  the  most 
basic  problems  in  graph  theory.  It  asks  for  the  shortest  path 
between  two  vertices  or  from  a  source  vertex  to  all  the 
other  vertices  (i.e.,  the  single-source  version  or  the  shortest 
path  tree).  Depending  on  whether  the  edge  weights  can  be 
negative,  the  problem  can  be  solved  via  Dijkstra’s  algorithm  or 
Bellman-Ford  algorithm  [2].  This  basic  problem  finds  diverse 
applications  in  communication  networks,  operational  research, 
plant  and  facility  layout,  and  VLSI  design  [3]. 

The  dynamic  version  of  the  shortest  path  problem  is  to 
maintain  the  shortest  path  tree  without  recomputing  from 
scratch  during  a  sequence  of  changes  to  the  graph.  A  typical 
change  to  a  graph  includes  weight  increase,  weight  decrease, 
edge  insertion,  and  edge  deletion.  The  last  two  types  of 
changes  model  network  topological  changes,  but  they  can  be 
conceptually  considered  as  special  cases  of  weight  changes  by 
allowing  weight  to  be  infinity.  Thus,  if  the  sequence  of  changes 
contains  only  weight  increase  and  edge  deletion,  we  call  it  a 
decremental  problem;  if  it  contains  only  weight  decrease  and 
edge  insertion,  we  call  it  an  incremental  problem.  Otherwise, 
we  have  a  fully  dynamic  problem.  If  multiple  edges  change 
simultaneously,  then  it  is  called  a  batch  problem. 

There  have  been  a  number  of  studies  of  the  dynamic 
shortest  path  problem  in  graphs.  Ramalingam  and  Reps  [4], 
Frigioni  et  al.  [5,  6],  and  Narvaez  et  al.  [7]  proposed 
several  algorithms  for  the  single-change  problem.  The  batch 
problem  was  considered  in  [7-9].  Comprehensive  experiments 
on  the  comparison  of  different  batch  algorithms  can  be  found 
in  [9,  10]. 

B.  Shortest  Path  Problem  in  Hypergraphs 

Both  the  static  and  dynamic  shortest  path  problems  have 
a  corresponding  version  in  hypergraphs.  The  static  shortest 
hyperpath  problem  was  considered  by  Knuth  [11]  and  Gallo 
et  al.  [12],  in  which  Dijkstra’s  algorithm  was  extended  to 
obtain  the  shortest  hyperpaths.  Knuth’s  algorithm  is  for  a 
special  class  of  hypergraphs  while  Gallo’s  algorithm  is  for 
a  general  hypergraph.  Ausiello  et  al.  proposed  a  dynamic 
shortest  hyperpath  algorithm  for  directed  hypergraphs,  con¬ 
sidering  only  the  incremental  problem  with  the  weights  of 
all  hyperedges  limited  to  a  finite  set  of  numbers  [13,  14],  A 
dynamic  algorithm  for  the  batch  problem  in  a  special  class  of 
hypergraphs  was  developed  in  [8]. 
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With  the  exception  of  the  above  few  studies,  the  shortest 
hyperpath  problem  remains  largely  unexplored.  To  the  best 
of  our  knowledge,  no  algorithms  exist  for  the  fully  dynamic 
problem  in  a  general  hypergraph. 

In  this  paper,  we  develop  two  fully  dynamic  shortest  path 
algorithms  for  general  hypergraphs.  These  two  algorithms 
complement  each  other,  with  each  preferred  in  different  types 
of  hypergraphs  and  dynamics. 

Referred  to  as  the  HyperEdge  based  Dynamic  Shortest  Path 
algorithm  (HE-DSP),  the  first  algorithm  is  an  extension  of  the 
dynamic  Dijkstra’s  algorithm  for  graphs  to  hypergraphs  (par¬ 
allel  to  Gallo’s  extension  of  the  static  Dijkstra’s  algorithm  to 
hypergraphs  in  [12]).  The  extension  of  the  dynamic  Dijkstra’s 
algorithm  to  hypergraphs  is  more  involved  than  that  of  the 
static  Dijkstra’s  algorithm.  This  is  due  to  the  loss  of  the  tree 
structure  (in  the  original  graph  sense)  in  the  collection  of  the 
shortest  hyperpaths  from  a  source  to  all  other  vertices.  Since 
the  dynamic  Dijkstra’s  algorithm  relies  on  the  tree  structure  to 
update  the  shortest  paths  after  an  incremental  change  (weight 
increase  or  edge  deletion),  special  care  needs  to  be  given  when 
extending  it  to  hypergraphs. 

The  second  algorithm  is  rooted  in  the  idea  of  Dimension 
Reduction  and  is  referred  to  as  DR-DSP.  The  basic  idea  is 
to  reduce  the  problem  to  finding  the  shortest  path  in  the 
underlying  graph  of  the  hypergraph.  The  underlying  graph  of  a 
hypergraph  has  the  same  vertex  set  and  has  an  edge  between 
two  vertices  if  and  only  if  there  is  at  least  one  hyperedge 
containing  these  two  vertices  in  the  original  hypergraph.  The 
weight  of  an  edge  in  the  underlying  graph  is  defined  as  the 
minimum  weight  among  all  hyperedges  containing  the  two 
vertices  of  this  edge.  The  shortest  hyperpath  in  the  hypergraph 
can  thus  be  obtained  from  the  shortest  path  in  the  underlying 
graph  by  substituting  each  edge  along  the  shortest  path  with 
the  hyperedge  that  lent  its  weight  to  this  edge.  The  correctness 
and  advantage  of  this  algorithm  are  readily  seen:  the  definition 
of  weight  in  the  underlying  graph  captures  the  minimum  cost 
offered  by  all  hyperedges  in  choosing  a  path  between  two 
vertices,  thus  ensuring  the  correctness  of  the  algorithm;  the 
reduction  of  a  hypergraph  to  its  underlying  graph  removes 
many  hyperedges  from  consideration  when  finding  the  shortest 
path,  leading  to  efficiency  and  agility  to  dynamic  changes. 

HE-DSP  is  more  efficient  in  hypergraphs  that  are  densely 
connected  through  high-dimensional  hyperedges  and  for  net¬ 
work  dynamics  where  changes  often  occur  to  hyperedges  that 
are  not  on  the  current  shortest  hyperpaths.  DR-DSP  has  lower 
complexity  when  hyperedge  changes  often  lead  to  changes  in 
the  shortest  hyperpaths.  This  is  usually  the  case  in  networks 
where  hyperedges  in  the  shortest  hyperpaths  are  more  prone 
to  changes  due  to  attacks,  frequent  use,  or  higher  priority  in 
maintenance  and  upgrade.  Furthermore,  DR-DSP  leads  to  an 
alternative  algorithm  for  solving  the  static  shortest  hyperpath 
problem  when  the  dynamic  problem  degenerates  to  the  static 
problem.  It  has  the  same  complexity  as  Gallo’s  algorithm 
for  a  general  hypergraph  and  lower  complexity  for  simplicial 
complexes  (a  special  class  of  hypergraphs  whose  hyperedge 
set  is  closed  under  the  subset  operation).  We  also  point  out 


that  both  proposed  algorithms  apply  to  directed  hypergraphs 
with  minor  modifications  in  their  implementation  details. 

A  detailed  time  complexity  analysis  of  these  two  algorithms 
is  provided  to  demonstrate  their  performance  in  the  worst- 
case  change  scenario.  Using  a  random  geometric  hypergraph 
model  and  a  real  data  set  of  a  social  network  (Enron  email 
data  set),  we  study  the  average  performance  of  these  two 
algorithms  in  different  scenarios  and  demonstrate  the  partition 
of  the  application  space  between  these  two  algorithms.  In  the 
experiment  with  Enron  email  data  set,  the  proposed  algorithms 
successfully  identified  the  most  important  actor  in  this  social 
network  using  the  closeness  centrality  metric. 

C.  Applications 

Shortest  path  computations  on  hypergraphs  can  be  applied 
to  communication  as  well  as  social  networks.  An  example 
application  in  wireless  communications,  in  particular,  for  mul¬ 
tihop  wireless  networks,  is  in  opportunistic  routing  schemes 
such  as  ExOR  [15],  GeRaF  [16],  and  MORE  [17].  In  such 
schemes,  any  receiver  of  a  packet  is  eligible  to  forward 
the  packet.  Receivers  typically  execute  a  protocol  amongst 
themselves  to  decide  who  should  forward  it.  This  naturally 
leads  to  a  hypergraph  model  where  a  node  and  its  neighbors 
form  a  hyperedge.  The  cost  of  each  hyperedge  can  be  defined 
based  on  the  cardinality  of  the  hyperedge  to  capture  the 
success  rate  of  forwarding  (lower  the  cardinality,  lesser  the 
chance  that  at  least  one  of  the  nodes  successfully  receives  the 
packet)  and  the  associated  overhead  (higher  the  cardinality, 
higher  the  energy  consumption  and  the  overhead  in  choosing 
the  forwarding  node).  A  shortest  hyperpath  from  the  source  to 
the  destination  is  thus  a  better  route  than  merely  the  traditional 
shortest  path.  And  as  the  network  topology  changes,  a  dynamic 
algorithm  is  required  to  maintain  the  shortest  hyperpath. 

In  social  networks,  information  (results,  event  reports,  opin¬ 
ions,  rumors,  etc.)  propagates  through  diverse  communication 
means  including  direct  links  (e.g.,  gestures,  optical,  satcom, 
regular  phone  call),  social  media  (e.g.,  Facebook,  Twitter, 
blogs),  mailing  lists,  and  newsgroups.  Such  a  network  may 
be  modeled  as  a  hypergraph  with  the  weight  of  a  hyperedge 
reflecting  the  cost,  credibility,  and/or  delay  for  disseminating 
information  among  all  vertice  of  this  hyperedge.  In  particular, 
the  weight  of  a  hyperedge  can  capture  the  unique  effect  on 
the  information  after  it  passes  through  a  group  of  people. 
For  instance,  a  result  can  be  discussed  by  overlapping  blog 
collaboration  networks  as  it  spreads,  and  often  the  discussion 
yields  a  better  result  than  if  it  only  spreads  through  individuals. 
The  minimum  cost  information  passing  in  social  networks  can 
thus  be  modelled  as  a  shortest  hyperpath  problem. 

Another  potential  application  is  that  of  finding  the  most 
important  actor  in  a  social  network.  Under  a  graph  model  of 
social  networks,  the  relative  importance  of  a  vertex  can  be 
measured  by  its  betweenness  and  closeness  centrality  indices. 
The  former  is  defined  based  on  the  number  of  shortest  paths 
that  pass  through  this  vertex,  and  the  latter,  the  total  weight  of 
the  shortest  paths  from  this  vertex  to  all  the  other  vertices  [18]. 
In  a  social  network  exhibiting  hyper-relationships,  between- 
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ness  and  closeness  centrality,  based  on  the  shortest  hyperpaths, 
would  be  better  indicators  of  the  relative  importance  of  each 
actor.  In  Sec.  VII,  we  apply  the  proposed  shortest  hyperpath 
algorithms  to  the  Enron  email  data  set.  We  propose  a  weight 
function  that  leads  to  the  successful  identification  of  the  CEO 
of  Enron  as  the  most  important  actor  under  the  closeness 
centrality  metric.  The  distance  of  each  person  in  the  data  set 
to  the  CEO  along  the  resulting  shortest  hyperpaths  closely 
reflects  the  position  of  the  person  within  the  company. 

II.  Background  on  Dynamic  Shortest  Path 
Algorithms  for  Graphs 

In  this  section,  we  present  the  basic  ideas  of  the  dynamic 
shortest  path  algorithms  developed  for  graphs  in  [5].  Some 
basic  techniques  in  updating  and  maintaining  the  shortest 
path  tree  will  be  borrowed  in  later  sections  when  we  develop 
dynamic  shortest  hyperpath  algorithms. 

A.  Dynamic  Shortest  Path  Problem 

A  change  6  on  a  graph  G  =  (V.  E )  corresponds  to  one 
edge  modification.  There  are  four  types  of  changes:  weight 
increase,  weight  decrease,  edge  insertion,  and  edge  deletion. 
Weight  increase  and  edge  deletion  can  be  similarly  treated 
(with  small  differences  in  the  required  data  structures  which 
will  be  omitted  for  simplicity),  so  can  weight  decrease  and 
edge  insertion.  The  dynamic  algorithms  are  thus  presented 
only  for  weight  increase  and  weight  decrease. 

Given  a  graph  G,  a  source  node  s,  and  a  sequence  of 
changes  C  =  {hi,  ^2,  •  •  • ,  Si}  on  G,  the  dynamic  shortest  path 
problem  is  to  find  the  shortest  paths  from  s  to  all  nodes  in 
each  new  graph  after  change  6, . 

In  the  following,  D[v]  denotes  the  distance  of  a  vertex  v  to 
the  source  s,  P\v\  the  parent  of  v  in  the  shortest  path  tree.  A 
vertex  v  is  called  an  affected  vertex  if  D[v\  or  P[v\  or  both 
change  in  the  new  shortest  path  tree.  An  edge  is  called  an 
affected  edge  if  it  contains  an  affected  vertex. 

B.  Weight  Decrease 

Consider  that  the  weight  of  edge  ( u,v )  decreases  to  wnew. 
Without  loss  of  generality,  assume  that  D[u\  <  D[v\.  It  is  not 
difficult  to  see  that  u  will  not  be  affected  by  this  change.  The 
dynamic  algorithm  starts  with  determining  whether  v  will  be 
affected  by  simply  checking  the  inequality 

D[u]  +  wnew  <  D[v}.  (1) 

If  the  inequality  does  not  hold,  then  this  edge  with  the 
decreased  weight  does  not  provide  a  shorter  path  for  v\  the 
algorithm  ends  and  the  shortest  path  tree  remains  unchanged. 
If  the  inequality  holds,  then  v  is  affected;  its  new  shortest 
path  from  s  must  go  through  edge  (u,v)  and  D[v ]  reduces  to 


D\u ]  +  wnew.  We  put  v  in  a  priority  queue1  Q,  and  the  rest 
of  the  procedure  is  similar  to  Dijkstra’s  algorithm:  dequeue 
the  node  2  with  the  minimum  distance  from  Q,  update  the 
distances  of  its  neighbors,  update  Q  by  inserting  the  new 
affected  vertices  among  the  neighbors  to  Q  and  update  the 
ranks  of  others  based  on  the  updated  distances.  The  procedure 
iterates  until  Q  is  empty.  A  pseudo  code  presentation  of  the 
basic  steps  is  given  below. 

Graph:  Weight  Decrease^,  v,  wnew). 

StepO  (Update  the  graph) 

1  w(u ,  V )  «-  Wnew 

Stepl  (Determine  the  affected  vertex  in  (u.  v'j) 

2  x<-  argmin qe{%i)}{D[q]};  y  «-  argma ^{aM{D\q\} 

3  if  D[x\  +  tV new  <  D[y)  do 

4  D[y]  <-  D[x\  +  wnew;  P[y\  <-  x 

5  Enqueue (Q,  (y,D[y])) 

6  end 

Step2  (Iteratively  update  all  affected  vertices) 

7  while  NonEmpty  (() )  do 

8  (z,D[z])  <r-  Dequeue(Q) 

9  for  each  v  £  V  s.t.  (z,v)  £  E 

10  if  D[v\  >  D[z]  +  w(z,v)  then 

11  D[v]  <—  D(z)  +  w(z,  v)\  P[u]  •£-  z 

12  Enqueue  or  Update  (Q,  (v,  D  [u])) 

13  end;  end;  end 

C.  Weight  Increase 

Consider  that  the  weight  of  edge  (u,v)  increases  to  wnew. 
Again,  assume  that  D[u\  <  D[v\.  If  (u,v)  is  not  an  edge 
in  the  shortest  path  tree,  then  none  of  the  vertices  will  be 
affected,  the  shortest  path  tree  remain  unchanged.  Otherwise, 
the  descendants,  and  only  the  descendants  of  this  edge  in  the 
shortest  path  tree  may  be  affected.  For  these  vertices,  some 
of  them  will  have  increased  distances,  some  of  them  will 
go  through  an  alternative  path  with  the  same  distance  (but 
changed  parent),  while  the  rest  will  not  be  affected.  In  order 
to  classify  the  vertices  into  these  three  categories,  we  introduce 
the  coloring  idea  in  Frigioni’s  algorithm  [5]: 

(1)  v  is  colored  white  if  neither  D[v\  nor  P[v\  needs  to  be 
changed. 

(2)  v  is  colored  pink  if  P[v]  needs  to  be  changed  but  D[v\ 
remains  the  same. 

(3)  v  is  colored  red  if  D[v\  increases. 

It  is  not  difficult  to  see  that  if  a  vertex  v  is  white  or  pink, 
all  its  descendants  in  the  shortest  path  tree  are  white;  if  v 
is  red,  all  its  descendants  are  either  red  or  pink.  Therefore 
the  coloring  procedure  is  clear:  we  first  determine  whether 
v  is  pink  or  red  by  checking  whether  there  is  an  alternative 

1 A  priority  queue  is  an  abstract  data  type  with  the  following  access  protocol: 
only  the  highest-priority  element  can  be  accessed.  Basic  operations  of  a 
priority  queue  include  Enqueue  (add  a  new  item  to  the  queue).  Dequeue 
(remove  the  item  with  the  highest  priority  and  return  this  item).  Update 
(change  the  priority  of  one  item  in  the  queue),  and  Peek  (obtain  the  value 
of  the  item  with  the  highest  priority).  Standard  implementations  of  a  priority 
queue  with  different  time  complexities  include  array,  link  list.  Binary  heap, 
and  Fibonacci  heap  [19]. 
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shortest  path  with  the  same  distance  for  v  (note  that  v  cannot 
be  white  due  to  the  weight  change  of  edge  (u,  v )  that  is  on 
its  current  shortest  path);  if  such  a  path  exists,  then  we  color 
v  pink  and  the  algorithm  ends,  otherwise  we  color  it  red  and 
put  all  its  children  in  a  priority  queue  M.  The  procedure  then 
iterates  for  each  vertex  in  M  according  to  an  increasing  order 
of  the  vertex  distances. 

After  the  coloring  process,  we  only  need  to  deal  with  the  red 
vertices.  For  each  red  vertex  z,  we  initialize  its  distance  with 
the  distance  of  the  shortest  path  through  one  of  its  non-red 
neighbors  and  put  z  in  another  priority  queue  Q  (if  no  non- 
red  neighbor  exists,  we  initialize  it  with  oo).  After  this,  the 
procedure  is  similar  to  Step  2  in  the  Graph:  Weight  Decrease 
algorithm:  at  each  iteration,  we  extract  the  vertex  at  the  top 
of  Q  and  update  its  neighbors  and  Q  until  Q  is  empty. 

Graph:  Weight  Increase^,  v,  wnew). 

StepO  (Update  the  graph) 

1  w{u,v)-^wnew 

Stepl  (Determine  the  affected  vertex  in  (u.  v )) 

2  x  <-  avgmmq&{ilM{D[q}} 

3  2/<-argmaxggKfi}{I%]} 

4  if  P[y]  =  x  then 

5  Enqueue(M,  (y,  D[y])) 

Step  2  (Coloring  Process) 

6  while  NonEmpty(M) 

7  (z,  D[z\)  -(—  Dequeue(M) 

8  if  3  nonred  q  £  V  s.t.  D[q]  +  w(q,  z)  =  D[z\ 

9  then  z  is  pink 

10  else  z  is  red;  Enqueue(M,  all  z’s  children) 

11  end; end 

Step3.a  (Initialize  the  distance  vector  for  red  vertices) 

12  for  each  red  vertex  z  do 

13  if  z  has  no  nonred  neighbor 

14  then  D[z\  < — boo;  P[z]  4—  Null 

15  else 

16  let  u  be  the  best  nonred  neighbor  of  z 

17  D[z\  4—  D[u]  +  w(u ,  z );  P[z]  4—  u 

18  Enqueue(Q,  (z,  D[z])) 

19  end; end;  end 

Step3.b:  Step2  of  Graph:  Weight  Decrease 

The  worst-case  time  complexity  for  one  edge  change  (either 
weight  decreasing  or  increasing)  is  0(|<5|  log  <5|  +  ||<5||),  where 
|b|  denotes  the  number  of  affected  vertices  and  ||5||  the  total 
number  of  both  affected  vertices  and  affected  edges. 

III.  Dynamic  Shortest  Hyperpath  Problem 

We  introduce  some  basic  concepts  of  hypergraph  [1]  and 
define  the  static  and  the  dynamic  shortest  hyperpath  problems. 
Some  basic  properties  of  the  shortest  hyperpaths  are  estab¬ 
lished  and  will  be  used  in  developing  the  dynamic  algorithms 
in  subsequent  sections. 


A.  Hypergraph  and  Hyperpath 

Let  The  a  finite  set  and  E  a  family  of  subsets  of  V.  If  for 
all  elements  et  £  E,  the  following  conditions  are  satisfied: 

G  7^  0;  ^-ieiGE  G  —  ^ 

then  the  couple  H  =  ( V.  E)  is  called  a  (undirected)  hyper¬ 
graph.  Each  element  v  €  V  is  called  a  vertex  and  each  element 
e  £  E  a  hyperedge. 

A  weighted  undirected  hypergraph  is  a  triple  H  =  ( V ,  E,  w) 
with  w  :  E  — {i?+  U  {0}}  being  a  nonnegative  weight 
function  defined  for  each  hyperedge  in  E. 

In  a  hypergraph,  a  hyperpath  is  defined  as  follows. 

Definition  1:  A  hyperpath  between  two  vertices  u  and  v  is 
a  sequence  of  hyperedges  {eo,  ei, . . . ,  em }  such  that  u  £  eo, 
v  £  em,  and  (T  g+i  ^  0  for  *  =  0, ...,  m  —  1.  A  hyperpath 
is  simple  if  non-adjacent  hyperedges  in  the  path  are  non¬ 
overlapping,  i.e.,  ei  D  ej  =  0,  Vj  ^  i,  i  ±  1. 

Let  Le  =  {eo,...,em}  be  a  hyperpath  in  a  weighted 
hypergraph  H.  We  define  the  weight  of  Le  as: 

m 

w(Le)  =  Ywjej). 

2=0 

B.  Shortest  Hyperpath  and  Relationship  Tree 

Given  two  vertices  u  and  v,  a  natural  question  is  to  find 
the  shortest  hyperpath  (in  terms  of  the  path  weight)  from  u 
to  v.  Since  the  weight  function  is  nonnegative,  it  suffices  to 
consider  only  simple  hyperpaths.  If  the  shortest  hyperpath  is 
not  simple,  we  can  always  generate  a  simple  hyperpath  without 
increasing  the  weight  by  deleting  all  the  hyperedges  between 
two  overlapping  non-adjacent  hyperedges. 

The  dynamic  shortest  hyperpath  problem  can  be  similarly 
defined  for  a  sequence  C  =  {bi,  <52,  •  •  ■ ,  of  hyperedge 
changes.  Hyperedge  changes  have  the  same  four  types  as 
edge  changes  in  a  graph:  weight  increase,  weight  decrease, 
hyperedge  insertion,  and  hyperedge  deletion.  Similarly,  weight 
increase  and  hyperedge  deletion  will  be  treated  together,  so  are 
weight  decrease  and  hyperedge  insertion. 

In  this  paper,  we  consider  the  single-source  shortest  hy¬ 
perpath  problem:  find  the  shortest  hyperpaths  from  a  given 
source  s  to  all  other  vertices.  The  presentation  of  the  paper 
focuses  on  undirected  hypergraphs.  However,  the  two  proposed 
dynamic  algorithms  apply  to  directed  hypergraphs  with  minor 
modifications  in  their  implementation  details. 

Below,  we  establish  a  basic  property  of  shortest  hyperpaths. 

Lemma  1:  Let  L  =  {e±,  e2,  . . . ,  e{\  be  a  shortest  hyperpath 
from  s  £  ei  to  z  £  e;.  Then  for  any  vertex  v  £  e*  D  e,+i,  the 
hyperpath  Lv  =  {ei,  e-j, ...  ■  fc, }  is  a  shortest  hyperpath  from 
s  to  v.  Furthermore,  for  any  two  vertices  u,  v  £  ei  FI  e^+i  (if 
there  exist  at  least  two  vertices  in  e*  n  e^+i),  D[u\  =  D[v\. 

Proof:  We  will  prove  by  contradiction.  Assume  that 
Lv  =  {ei,  e2,  •  ■  • ,  ei}  is  not  a  shortest  hyperpath  for  v. 
Then  there  exists  a  different  hyperpath  If,  =  {e\ ,  e'2. ....  e'k } 
with  w(L'v )  <  w(Lv).  Then  consider  the  hyperpath  L'  = 
{ei,e'2,...,e'fe,ei+i ,  ei+2, .. .  ,e*},  we  have  w(L')  <  w{L) 
which  contradicts  the  fact  that  L  is  a  shortest  hyperpath  to 
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z.  This  completes  the  proof  for  the  first  part  of  the  lemma. 
Furthermore,  for  any  two  nodes  it,  v  &  e^nei+i,  since  Lv  is  the 
shortest  hyperpath  for  both  vertices,  D[v\  =  w(Lv)  =  D[u ]. 

■ 

Next,  we  introduce  the  concept  of  relationship  tree  that  is 
needed  in  the  proposed  dynamic  shortest  hyperpath  algorithm 
HE-DSP.  Since  two  adjacent  hyperedges  in  a  hyperpath  may 
overlap  at  more  than  one  vertex,  the  shortest  hyperpaths  from 
s  to  all  other  vertices  do  not  generally  form  a  tree  in  the 
original  graph  sense.  For  the  development  of  the  dynamic 
shortest  hyperpath  algorithms,  we  introduce  the  concept  of 
relationship  tree  to  indicate  the  parent-child  relationship  along 
shortest  hyperpaths.  The  concept  can  be  easily  explained  in  the 
example  given  in  Fig  2.  Let  {ei,e2}  be  a  shortest  hyperpath 
from  s  to  ('4 .  By  Lemma  1,  {e±}  is  a  shortest  hyperpath  for 
both  v\  and  v-2-  As  illustrated  in  Fig  2,  there  are  4  possible 
relationship  trees  to  indicate  the  parent-child  relationship  in 
these  shortest  hyperpaths.  We  will  show  in  Sec.  IV  that  the 
choice  of  the  relationship  tree  does  not  affect  the  correctness 
or  performance  of  the  proposed  algorithm  HE-DSP. 


Fig.  2.  Hyperpaths  and  the  associated  relationship  trees. 


Similar  notations  are  used  for  dynamic  shortest  hyperpath 
algorithms:  D[v ]  denotes  the  distance  of  a  vertex  v  to  the 
source  s  on  the  shortest  hyperpath,  P[v]  the  parent  of  v 
in  the  chosen  relationship  tree  associated  with  the  shortest 
hyperpaths.  A  new  notation  is  E[v],  the  hyperedge  containing 
v  and  P[v\  on  the  shortest  hyperpath  (i.e.,  the  hyperedge  that 
leads  to  v  from  P\v]  on  the  shortest  hyperpath).  When  it  is 
necessary  to  distinguish  the  shortest  distance  before  and  after 
a  weight  change,  d[v]  denotes  the  shortest  distance  before 
the  change,  d![v\  the  shortest  distance  after  the  change,  and 
D[v\  the  actual  value  stored  in  the  data  structure  during  the 
execution  of  the  algorithm. 

IV.  Hyperedge  Based  Dynamic  Shortest  Path 
Algorithm 

In  this  section,  we  propose  HE-DSP.  It  is  an  extension  of  the 
dynamic  Dijkstra’s  algorithm  to  hypergraphs.  The  extension  is 
more  complex  than  Gallo’s  extension  of  the  static  Dijkstra’s 
algorithm,  since  the  dynamic  Dijkstra’s  algorithm  relies  on  the 
tree  structure  of  the  shortest  paths,  a  structure  no  longer  there 
for  the  shortest  hyperpaths. 


A.  Hyperedge  Weight  Decrease 

Consider  that  the  weight  of  a  hyperedge  e  decreases  to 
wnew.  Similar  to  the  case  for  graphs,  we  know  that  the  vertex 
x  £  e  with  D[ x\  =  min„e,s {£)[?;]}  will  not  be  affected.  We 
then  check  weather  the  other  vertices  in  e  are  affected  by 
checking  the  inequality  given  in  (1),  and  put  all  the  affected 
vertices  into  a  priority  queue  Q.  The  rest  of  the  procedure  is 
similar  to  that  for  graphs,  only  when  we  update  the  distance  of 
a  vertex,  we  check  all  the  hyperedges  that  contain  this  vertex. 

HE-DSP:  Weight  Decrease  (  e,  wnew). 

StepO  (Update  the  hypergraph) 

1  w(e )  <-  wnew 

Stepl  (Determine  the  affected  vertices  in  e) 

2  x<<-argmin  vee{D[v}} 

3  for  each  v  G  e  such  that  D[x]  +  wnew  <  D[v]  do 

4  D[v\  <—  D[ x\  +  wnew\  P[v\  t—  x;  E[v\  <—  e 

5  Enqueue(Q,  ( v ,  D[v])) 

6  end 

Step2  (Iteratively  enqueue  and  update  affected  vertices) 

7  while  NonEmpty  (Q)  do 

8  (z,  D[z])  <—  Dequeue(Q) 

9  for  each  e  €  E  s.t.  z  £  e 

10  for  each  v  G  e 

11  if  D[v\  >  D[z\  +  w(e)  then 

12  D[v\  <—  D(z)  +  w(e);  P[v ]  z;  E[v]  -(—  e 

13  Enqueue  or  Update(Q,  (v,  -D[i>])) 

14  end;  end;  end;  end 

Theorem  1:  If  before  the  weight  decrease,  D[v\  =  d[r>], 
E[v]  and  P[v]  are  correct  for  all  v  £  V,  then  after  the 
weight  decrease,  D[v\  =  d'[v\  and  E[v\  and  P\v\  are  correctly 
updated. 

Proof:  See  Appendix  A.  ■ 

B.  Hyperedge  Weight  Increase 

The  coloring  process  in  the  graph  case  relies  on  the  tree 
structure  of  the  shortest  paths,  which  is  no  longer  present  in 
the  shortest  hyperpaths.  Our  solution  is  to  use  a  relationship 
tree  for  the  coloring  process,  and  we  prove  the  correctness  of 
this  approach  regardless  of  the  choice  of  the  relationship  tree. 

Consider  that  the  weight  of  a  hyperedge  e  increases  to  wnew . 
First,  we  redefine  the  color  of  a  vertex  v  based  on  the  chosen 
relationship  tree. 

(1)  v  is  colored  white  if  d'[v\  =  d[v\  while  keeping  the 
current  P[v]  and  E[v\. 

(2)  v  is  colored  pink  if  d'[v]  =  d[r>],  but  only  possible  through 
a  new  P[v]  or  E[v)  or  both. 

(3)  v  is  colored  red  if  d'[v\  <  d[v\. 

With  the  above  modified  definitions  of  colors,  the  same 
coloring  process  as  in  the  graph  case  can  be  carried  out  using 
a  relationship  tree.  The  algorithm  is  given  below. 

HE-DSP:  Weight  Increase(e,  wnew). 

StepO  (Update  the  hypergraph) 

1  w(e)  <-  wnew 

Stepl  (Determine  the  affected  vertices  in  e) 
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2  for  each  v  G  e  s.t.  E[v]  =  e  do 

3  Enqueue(M,  (v,  D[v})) 

Step2  (Coloring  process) 

4  while  NonEmpty(M) 

5  (z,  D[z])  G-  Dequeue(M) 

6  if  3  nonred  q  £  V  s.t.  3 e  G  E  with  q,z  G  e  and 

D[q\  +  w(e)  =  D[z ) 

7  then  z  is  pink;  P[z]  =  q\  E[z\  =  e; 

8  else  z  is  red;  Enqueue(M,  all  z’s  children) 

9  end;  end 

Step3.a  (Initialize  the  distance  vector  for  red  vertices) 

10  for  each  red  vertex  z  do 

11  if  z  has  no  nonred  neighbor 

12  then  D[z]  < — boo;  P[z]  G-  Null 

13  else 

14  let  it  be  the  best  nonred  neighbor  of  z 

15  E[z\  G-  argminegB  e3u  ,{w(e)}; 

16  D[z\  G-  D[u ]  +  w(E[z\);  P[z ]  G-  u; 

1 7  Enqueue (Q ,  (z,D[z})) 

18  end; end; end 

Step3.b:  Step2  of  HE-DSP:  Weight  Decrease 

The  theorem  below  states  the  correctness  of  the  algorithm. 
Theorem  2:  If  before  the  weight  increase,  D[v\  =  <7[ii], 
E[v]  and  P[t>]  are  correct  for  all  v  G  V,  then  after  the  weight 
increase,  D[v\  =  d'[v\  and  also  E[v\  and  P[t;]  are  correctly 
updated. 

Proof:  See  Appendix  B.  ■ 

V.  Dimension  Reduction  based  Dynamic  Shortest 
Path  Algorithm 

In  this  section,  we  propose  DR-DSP.  When  the  dynamic 
problem  degenerates  to  the  static  problem,  DR-DSP  leads  to  an 
alternative  algorithm  for  solving  the  static  shortest  hyperpath 
problem. 

A.  The  Static  Case:  DR-SP 

We  first  consider  the  static  version  of  the  algorithm  (referred 
to  as  DR-SP),  which  captures  the  basic  idea  of  dimension 
reduction. 

The  proposed  DR-SP  algorithm  is  based  on  the  following 
theorem  in  which  we  show  that  for  a  general  hypergraph  H, 
the  weight  oj(L*)  of  the  shortest  path  L*  of  II  is  equal  to 
the  shortest  path  L*G  of  a  weighted  graph  G  derived  from 
H.  Specifically,  corresponding  to  every  hyperedge  e  in  77,  G 
contains  a  clique  defined  on  the  vertices  of  e. 

Theorem  3:  Let  77  =  (V,  E ,  w)  be  a  hypergraph,  and  G  = 
(V,  E)  the  underlying  graph  of  77  where  an  edge  e  G  73  if 
and  only  if  3e  G  E  such  that  e  C  e.  For  each  edge  e  in  G,  its 
weight  wc ?(e)  is  defined  as 

wo(e)  =  min  w(e).  (2) 

{eE.E:  eDe} 

Let  L*  and  LG  be  the  shortest  paths  from  u  G  V  to  v  G  V  in 
77  and  G,  respectively.  Then  we  have  that 

w(L*)  =  wg(L*g). 


Proof:  First,  for  each  shortest  path  L*G  in  G,  we  can 
obtain  a  corresponding  hyperpath  L  in  77  with  the  same  weight 
based  on  (2),  therefore  we  have  that 

wg(Lg)  =  w(L)  >  w(L*). 

Then  it  suffices  to  show  that  there  exists  a  path  Lq  in  G 
such  that  wg{Lq )  <  w{L*),  which  implies  that  wg{Lg)  < 
wg(Lg)  <  w(L*). 

Assume  that  L*  =  {eo,  e±, . . . ,  ek-ij  is  a  shortest  hyper¬ 
edge  path  from  vo  to  Vk  in  H  where  vq  G  eo  and  Vk  G  ek-i- 
Let  Vi  G  ej_ i  (~l  et  (i  =  1, 2, ...,  k  —  1)  be  one  of  the  vertices 
in  the  intersection  of  hyperedges  ej_i  and  e,.  Construct  a 
path  Lg  =  {vq,Vi ,...,Vk}  in  the  graph  G.  For  each  edge 
ei  =  {ui,  r’i+i}  (i  =  0, 1, ...,  k  —  1),  since  e,  C  a,  it  follows 
from  (2)  that 

wcief)  <  w{ef). 

Thus, 

k- 1  k- 1 

wg(Lg )  =  ^2  wg{g)  <  ^2  w (ei)  = 

7—0  7—0 

i.e.,  wq{Lg)  <  w{L*).  ■ 

It  follows  from  Theorem  3  that  the  shortest  path  in  a  general 
hypergraph  can  be  obtained  by  applying  Dijkstra’s  algorithm 
to  the  underlying  graph  G  with  weights  modified  as  stated  in 
the  theorem. 

B.  The  Dynamic  Case:  DR-DSP 

In  the  dynamic  case,  a  sequence  G  =  {<5i,  62,  ■  ■  ■ ,  Si}  of 
hyperedge  changes  in  the  hypergraph  77  results  in  a  sequence 
of  edge  changes  in  the  underlying  graph  G.  For  each  hyper¬ 
edge  change  <5?:,  DR-DSP  first  updates  the  underlying  graph 
G  to  locate  all  the  changed  edges  caused  by  5i.  In  the  next 
step,  DR-DSP  updates  the  shortest  path  tree  in  the  underlying 
graph  G. 

Consider  first  the  graph  update.  A  change  to  a  hyperedge 
e  only  affects  those  edges  in  G  that  are  subsets  of  e,  i.e., 
a  hyperedge  change  is  localized  in  the  underlying  graph  G. 
Furthermore,  since  the  weight  of  an  edge  in  G  is  the  minimum 
weight  of  all  hyperedges  containing  it,  not  all  edges  in  G  that 
are  subsets  of  e  will  change  weight.  Based  on  these  obser¬ 
vations,  we  propose  a  special  data  structure  and  procedure 
for  updating  the  underlying  graph  G  without  regenerating  the 
graph  from  scratch  using  Step  1  of  DR-SP. 

At  the  initialization  stage  of  the  algorithm,  a  priority  queue 
Muv  for  each  pair  of  vertices  (it,  v)  in  the  hypergraph  is 
established  to  store  the  weights  of  all  hyperedges  that  contain 
both  u  and  v.  When  a  change  occurs  to  hyperedge  e,  all  the 
priority  queues  Muv  associated  with  the  pair  of  vertices  (it,  v) 
that  are  contained  in  e  are  updated  with  the  new  weight  of  e. 
Thus,  the  top  of  these  priority  queues  always  maintain  the 
weight  for  edge  (it,  v)  in  the  underlying  graph  G  for  each 
(it,  v).  Below  is  a  pseudo  code  implementation  of  the  proposed 
procedure. 

Graph  Update(e,  wnew). 
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1  for  each  it,  v  €  e 

2  Updat e(Muv,  <  e,  wnew  >); 

3  wuv  <-Peek(M„w); 

4  end; 

After  the  underlying  graph  G  is  updated,  we  are  now  facing 
a  dynamic  shortest  path  problem  in  a  graph.  However,  since  a 
single  hyperedge  change  can  result  in  multiple  edge  changes 
in  G,  we  need  to  handle  a  batch  problem.  While  existing 
batch  algorithms  and  iterative  single-change  algorithms  for 
graphs  can  be  directly  applied  here,  we  show  that  the  batch 
problem  we  have  at  hand  has  two  unique  properties  that  can 
be  exploited  to  improve  the  efficiency  of  the  algorithm. 

Property  1:  The  edge  changes  in  G  caused  by  a  hyperedge 
change  are  either  all  weight  decreases  or  all  weight  increases. 

Property  2:  All  changed  edges  in  G  caused  by  a  hyperedge 
change  belong  to  a  clique  in  G. 

C.  Hyperedge  Weight  Decrease 

If  the  weight  of  hyperedge  e  decreases  to  wnew  ,  by 
Theorem  3  and  Property  1,  there  are  (possibly)  several  edge- 
weight  decreases  in  the  underlying  graph  G.  Therefore  sim¬ 
ilar  to  HE-DSP,  there  is  at  least  one  unaffected  node  x  = 
argminued{D[i;]}.  By  Property  2,  these  affected  edges  are 
contained  in  a  clique  derived  from  the  changed  hyperedge; 
therefore  it  is  sufficient  to  determine  the  distance  of  every 
node  v  (other  than  x)  in  the  original  changed  hyperedge  e 
by  checking  D[x\  +  wnew  <  D[v\.  And  we  can  initialize  the 
priority  queue  with  those  nodes  whose  weight  decreases.  After 
that,  the  procedure  is  similar  to  that  in  the  graph  case. 

DR-DSP:  Weight  Decrease  (  e,  wnew). 

StepO  (Update  the  hypergraph  and  G ) 

1  w(e)  <-  wnew 

2  Graph  Update(e,  wnew) 

Stepl  of  HE-DSP:  Weight  Decrease 

Step2  of  Graph:  Weight  Decrease 

D.  Hyperedge  Weight  Increase 

If  the  weight  of  hyperedge  e  increases  to  wnew,  by  Theo¬ 
rem  3  and  Property  1,  there  are  (possibly)  several  edge-weight 
increases  in  the  underlying  graph  G.  Similar  to  the  single¬ 
change  case  in  graph,  there  is  at  least  one  unaffected  node 
x  =  argmin7,Cf, {/J[u] }.  Then  another  node  v  £  e  is  affected 
only  if  E[v\  =  e,  i.e.,  e  is  on  its  shortest  hyperpath.  We  use 
all  such  nodes  to  initialize  the  priority  queue  M.  The  rest  is 
similar  to  the  procedure  of  Graph:  Weight  Increase. 

DR-DSP:  Weight  Increase (  e,  wnew). 

StepO  (Update  the  hypergraph  and  G ) 

1  w(e)  <-  wnew 

2  Graph  Update(e,  wnew) 

Stepl  of  HE-DSP:  Weight  Increase 

Step2  of  Graph:  Weight  Increase 

Step3.a  of  Graph:  Weight  Increase 

Step3.b  of  Graph:  Weight  Increase 


VI.  Time  Complexity  Analysis 

We  analyze  the  time  complexity  of  the  two  proposed  dy¬ 
namic  algorithms.  We  show  that  for  different  scenarios,  each 
algorithm  has  its  own  advantage.  We  also  consider  the  static 
case  and  show  that  the  static  version  of  DR-DSP  has  the  same 
complexity  as  Gallo’s  algorithm  for  a  general  hypergraph  and 
lower  complexity  for  a  simplicial  complex. 

A.  The  Static  Shortest  Hyperpath  Problem 

Given  a  hypergraph  H  =  (V,  E,  w),  let  n  =  \V\  denote  the 
number  of  vertices  in  H,  and  $  =  YheeE  le|2  where  |e|  is 
the  cardinality  of  e.  For  a  simplicial  complex,  let  m  be  the 
number  of  facets,  and  d  the  maximum  degree  of  the  facets. 

Theorem  4:  The  time  complexities  of  Gallor’s  algorithm 
and  DR-SP  for  general  hypergraphs  and  simplicial  complexes 
are  as  follows. 


Algorithm 

General  Hypergraph 

Simplicial  Complex 

Gallo 

0(n  log  n  +  <£) 

0(n  log  n  +  d/!2dm) 

DR-SP 

0(n  log  n  +  <h) 

0(n  logn  +  d2dm ) 

Proof:  The  time  complexity  of  DR-SP  mainly  comes 
from  Steps  1  and  2.  Step  2  is  essentially  applying  Dijkstra’s 
algorithm  to  a  graph  with  n  vertices  and  rh  edges  where  rh  is 
the  number  of  edges  in  the  underlying  graph  G.  The  running 
time  is  thus  0(n\ogn  -pin).  An  implementation  of  Step  1  is 
to  obtain  the  edge  weight  wc{e)  based  on  (2).  Therefore  the 
time  complexity  for  Step  1  is  0(fffeeE  |e|2),  i.e.,  0(<t>).  With 
rh  upper  bounded  by  $  (since  for  each  e  €  E,  there  are  at 
most  |e|(|e|  —  l)/2  edges  in  G),  we  arrive  at  the  total  time 
complexity  of  DRSP 

For  Gallo’s  Algorithm,  similar  to  Dijkstra’s  algorithm,  the 
time  complexity  is  mainly  in  updating  the  neighbors  of  the 
non-fixed  vertex  z  with  the  minimal  distance  D\z\.  For  each 
z,  the  algorithm  scans  all  the  hyperedges  containing  z.  For 
each  pair  of  vertices  (it,  v)  £  e,  e  is  scanned  twice.  Therefore 
the  total  number  of  such  operations  is  $  =  XXeB  lel2- 
Also,  extracting  2  from  the  priority  queue  implemented  by  a 
fibonacci  heap  takes  O(logn)  time.  The  total  time  complexity 
of  Gallo’s  algorithm  thus  follows. 

For  a  simplicial  complex,  $  =  0(d22dm),  the  complexity 
of  Gallo’s  algorithm  thus  follows.  For  DR-SP,  exploiting  the 
property  that  the  edge  set  is  closed  under  the  subset  operation 
in  a  simplicial  complex,  we  can  use  a  top-down  scheme  in 
Step  1  of  DR-SP  to  calculate  the  weight  vjg(s )  inductively 
with  respect  to  the  dimension  of  a  facet  as  follows: 

w g{s)  =  min{u;(s),  {iug^OI  s'  D  s  and  dim[s']  =  i  +  1}}, 

where  wg{s')  =  w(s')  for  the  facet  s'.  The  time  complexity 
for  Step  1  can  then  be  improved.  Because  each  i— dimensional 
face  is  associated  with  d  —  i  comparisons.  Thus,  the  running 
time  of  Step  1  for  each  d-dimensional  facet  is  given  by 

d- 1 

i—  1 


Therefore  the  time  complexity  for  Step  1  is  0(d2dm).  The 
total  time  complexity  thus  follows.  ■ 

B.  The  Dynamic  Shortest  Hyperpath  Problem 

Given  a  hypergraph  H  =  (V.  E,  w )  and  a  change  to 
hyperedge  e,  let  |<S|  denote  the  number  of  affected  vertices, 

1 1 <5 1 1  the  number  of  affected  hyperedges  plus  |5|,  |<5$|  = 
EeeB.e  is  affected  |e|2’  and  ll^ll  the  number  of  affected  edges 
in  the  underlying  graph  plus  |5|. 

Theorem  5:  The  time  complexities  of  HE-DSP  and  DR- 
DSP  for  the  fully  dynamic  shortest  path  problem  in  a  general 
hyperpath  are  as  follows. 


Algorithm 

Time  Complexity 

HE-DSP 

o(Hiog|a|  +  |A-*|) 

DR-DSP 

0(|<5|  log  |*|  +  ||*||  +  \e\2  log m) 

Proof:  For  HE-DSP:  Weight  Decrease,  the  dominating 
part  is  Step  2.  In  Step  2,  there  are  total  |5|  iterations.  In  each 
iteration,  the  algorithm  first  dequeues  one  node  z  from  M 
which  takes  0(log  |6|)  time.  Then  the  algorithm  updates  all 
of  z’s  neighbors  by  scanning  all  the  hyperedges  containing  2. 
Each  affected  hyperedge  e  can  be  scanned  at  most  |e|(|e|  — 
1)  =  0(|e|2)  times.  Therefore  the  time  spent  on  updates  for  all 
iterations  is  0(|<5$|).  The  total  time  complexity  of  HE-DSP: 
Weight  Decrease  is  0(|<5|  log  |<5|  +  |5$|).  For  HE-DSP:  Weight 
Increase,  similar  to  the  above  analysis,  the  time  complexity 
for  the  dominating  part  (Step  2,  Step  3. a  and  Step  3.b)  is 
0(|i5|  log  |(5|  +  |<5$|).  The  total  time  complexity  of  HE-DSP: 
Weight  Increase  is  0(|<5|  log  |5|  +  |5$|).  The  result  thus  follows. 

For  DR-DSP:  Weight  Decrease,  the  total  time  spent  on 
Graph  Update  procedure  is  0(|e|2  log rn).  In  Step  2,  there 
are  |$|  iterations;  in  each  iteration  0(log  |<5|)  time  is  spent  to 
dequeue  z  from  M.  Time  spent  on  updating  neighbors  over 
all  iterations  is  0(||5||).  Therefore  the  total  time  complexity  is 
0(|<5|  log  |()|  +  ||(5||-|-|e|2  log  to).  For  DR-DSP:  Weight  Increase, 
Step  2,  Step  3. a  and  Step  3.b  take  0(|<5|  log  |<5|  +  ||<5||)  (similar 
to  the  analysis  for  graphs).  The  total  time  complexity  thus 
follows.  ■ 

From  Theorem  5  we  see  that  if  <5  is  small  and  |e|  is 
large,  HE-DSP  performs  better,  since  in  DR-DSP,  the  update 
of  the  underlying  graph  has  to  be  done  regardless  whether 
there  are  affected  vertices.  Thus  in  a  sequence  of  hyperedge 
changes,  if  only  a  small  fraction  of  them  actually  have  affected 
nodes,  then  HE-DSP  will  outperform  DR-DSP.  On  the  other 
hind,  if  6  is  large,  because  usually  |<5$|  ||<5||,  then  DR- 

DSP  will  outperform  HE-DSP.  Consider  the  extreme  example 
where  every  valid  hyperedge  exists,  all  nodes  are  affected  and 
the  changed  hyperedge  contains  n  vertices.  Then  |d  =  n, 
|fo|  =  0{n2 2"),  |5|  =  0(n2),  |e|  =  n,  m  =  0(2").  The 
time  complexity  of  HE-DSP  is  0(nlogn  +  n22")  =  0(n2 2") 
while  the  time  complexity  of  DR-DSP  is  0(n3).  We  see  that 
the  time  complexity  of  DR-DSP  can  be  much  lower  than  that 
of  HE-DSP. 


VII.  Simulation  Results 

We  present  simulation  results  on  the  running  time  of  the 
proposed  dynamic  shortest  hyperpath  algorithms.  We  test  the 
proposed  algorithms  on  hypergraphs  generated  from  a  random 
geometric  model  as  well  as  those  generated  by  the  Enron  email 
data  set.  All  simulation  code  is  compiled  and  run  on  the  same 
laptop  equipped  with  a  3.0GHz  i7-920XM  Mobile  Processor. 

A.  Random  Geometric  Hypergraph 

We  first  consider  a  random  geometric  hypergraph  model  in 
which  n  nodes  are  uniformly  distributed  in  an  a  x  a  square.  All 
nodes  within  a  circle  with  radius  r  form  a  hyperedge  (circles 
are  centered  on  a  h  x  h  grid).  The  weight  of  each  hyperedge 
is  given  by  the  average  distance  between  all  pairs  of  vertices 
of  this  hyperedge. 

A  sequence  of  changes  are  then  generated  and  the  proposed 
dynamic  algorithms  are  employed  to  maintain  all  the  shortest 
hyperpaths  from  the  source  s  located  at  a  corner  of  the 
ax  a  square.  Each  change  can  be  a  hyperedge  insertion  (with 
probability  pj),  a  hyperedge  deletion  (with  probability  Pd),  or 
a  weight  change  (with  probability  1  —pi—pif)  with  new  weight 
chosen  uniformly  in  [vjrnin .  wmax].  In  the  case  of  a  hyperedge 
deletion  or  a  weight  change,  the  hyperedge  to  be  deleted  or  to 
be  assigned  with  a  new  weight  is  chosen  according  to  the 
two  models  detailed  below.  Hyperedge  insertions  are  only 
realized  when  there  are  hyperedges  that  have  been  deleted, 
and  a  randomly  chosen  one  is  inserted  back.  This  ensures  that 
all  hyperedges  satisfy  the  geometric  property  determined  by 
r  at  all  time.  It  also  models  the  practical  scenario  where  a 
broken  link  is  repaired. 

In  selecting  a  hyperedge  for  deletion  or  weight  change, 
we  consider  two  different  models:  the  random  change  model 
and  the  targeted  change  model.  In  the  former,  the  hyperedge 
is  randomly  and  uniformly  chosen  among  all  hyperedges. 
In  the  latter,  it  is  randomly  and  uniformly  chosen  from  the 
current  shortest  hyperpaths.  This  models  the  scenarios  where 
hyperedges  in  the  shortest  hyperpaths  are  more  prone  to 
changes  due  to  attacks,  frequent  use,  or  higher  priority  in 
maintenance  and  upgrade. 

In  Fig.  3,  we  show  the  simulation  results  on  the  running 
time  of  the  two  proposed  algorithms  under  a  sequence  of 
104  changes.  We  see  that  HE-DSP  has  lower  complexity  in 
networks  with  random  topological  and  weight  changes  (Fig.  3- 
Left),  whereas  DR-DSP  should  be  preferred  in  networks 
with  targeted  changes  (Fig.  3-Right).  This  partition  of  the 
application  space  can  be  explained  from  the  structures  of  these 
two  algorithms.  Under  the  random  change  model,  a  large 
fraction  of  changes  do  not  result  in  changes  in  the  current 
shortest  hyperpaths.  Such  changes  lead  to  little  computation 
in  maintaining  the  shortest  hyperpaths  for  both  algorithms, 
but  requires  about  the  same  amount  of  computation  in  the 
Graph-Update  step  of  DR-DSP  for  maintaining  the  underlying 
graph.  On  the  other  hand,  under  the  targeted  change  model, 
all  hyperedge  deletions  and  weight  changes  affect  the  shortest 
hyperpaths.  Updating  the  shortest  hyperpaths  can  be  done 
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more  efficiently  in  DR-DSP  since  it  works  on  the  underlying 
graph  with  a  much  smaller  number  of  edges. 


Fig.  3.  The  average  running  time.  Left:  the  random  change  model;  Right: 
the  targeted  change  model  (n  =  1000,  a  =  1000,  r  =  VT000,  h  =  1, 
Pi  =  Pd  —  ^min  —  10,  tUmax  =  20,  the  average  is  taken  over  50 
random  hypergraphs). 


2.5 


President  VP  MD  Manager  Director  Trader  Employee 


Fig.  4.  The  average  distance  from  the  CEO  to  others  at  different  positions. 


B.  Enron  Email  Data  Set 

In  this  example,  we  consider  the  application  of  the  shortest 
hyperpath  algorithms  in  finding  the  most  important  actor  in  a 
social  network.  We  consider  the  Enron  email  data  set  and  use 
the  same  hypergraph  generation  model  as  in  [20].  Specifically, 
each  person  is  a  vertex  of  the  hypergraph,  and  the  sender  and 
recipients  of  every  email  form  a  hyperedge.  Our  objective  is  to 
identify  the  most  important  person  measured  by  the  closeness 
centrality  index  (i.e.,  the  total  weight  of  the  shortest  hyperpaths 
from  this  person  to  all  the  other  persons).  The  first  step  is  to 
assign  weight  to  each  hyperedge  that  reflects  “distance”.  While 
there  is  no  universally  accepted  way  of  measuring  distance  in 
a  social  network  observed  through  email  exchanges,  certain 
general  rules  apply.  First,  a  direct  email  exchange  between 
two  persons  indicates  a  stronger  tie  than  an  email  sent  to 
a  large  group.  Thus,  the  weight  of  an  hyperedge  should  be 
an  increasing  function  of  the  cardinality  of  this  hyperedge. 
Second,  more  frequent  email  exchange  among  a  given  group  of 
people  shows  stronger  ties.  Thus,  the  weight  of  an  hyperedge 
should  be  decreasing  with  the  number  of  times  that  this 
hyperedge  appears  in  the  email  data  set.  Considering  these 
two  general  rules,  we  adopt  the  following  weight  function: 

w(e)  =  (a/N)“(1  X)  (3) 

where  \e\  is  the  cardinality  of  the  hyperedge  e,  a  is  the 
parameter  measuring  how  fast  the  weight  decreases  with  the 
number  l  of  times  that  this  hyperedge  appears  in  the  data  set. 

We  can  then  apply  DR-SP  on  the  resulting  (static)  hyper¬ 
graph  to  find  the  shortest  hyperpaths  rooted  at  each  vertex 
and  compute  this  vertex’s  closeness  centrality  index.  With  the 
weight  function  given  in  (3)  using  a  =  0.6,  the  identified  most 
important  actor  is  the  CEO  of  Enron.  The  average  distance 
(along  the  shortest  hyperpath)  from  the  CEO  to  the  other 
persons  at  various  positions  is  shown  in  Fig.  4.  We  observe 
that  in  general,  the  higher  the  position,  the  shorter  the  distance. 
These  results  demonstrate  that  the  adopted  hypergraph  model 
and  weight  function  capture  the  essence  of  the  problem. 

Next,  we  construct  a  dynamic  hypergraph  sequence  based 
on  the  Enron  data  set.  At  the  beginning,  the  hypergraph 


contains  only  individual  vertices.  We  then  consider  each  email 
chronologically.  Each  email  either  adds  a  new  hyperedge  or 
decrease  the  weight  of  an  existing  hyperedge  (due  to  the 
increased  number  of  appearances  of  this  hyperedge).  The  two 
proposed  algorithms  are  employed  to  maintain  the  shortest 
hyperpaths  rooted  at  the  CEO  after  each  change.  The  running 
time  is  given  in  Fig.  5,  which  shows  the  lower  complexity 
of  DR-DSR  The  reason  is  that  a  large  fraction  of  hyperedge 
changes  result  in  changes  in  the  shortest  hyperpaths. 


number  of  hyperedge  changes  x  104 

Fig.  5.  The  average  running  time  for  the  Enron  data  set  (a  =  0.6,  the 
average  is  taken  over  50  monte  carlo  runs). 

VIII.  Conclusion 

We  have  presented,  to  our  best  knowledge,  the  first  study 
of  the  fully  dynamic  shortest  path  problem  in  a  general  hyper¬ 
graph.  We  have  developed  two  dynamic  algorithms  for  finding 
and  maintaining  the  shortest  hyperpaths.  These  two  algorithms 
complement  each  other  with  each  one  preferred  in  different 
types  of  hypergraphs  and  network  dynamics,  as  illustrated  in 
the  time  complexity  analysis  and  simulation  experiments.  We 
have  discussed  and  studied  via  experiments  over  a  real  data 
set  the  potential  applications  of  the  dynamic  shortest  hyperpath 
problem  in  social  and  communication  networks. 

Appendix  A:  Proof  of  Theorem  1 
The  proof  is  based  on  the  following  three  lemmas. 
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Lemma  2:  Let  x  =  argmin^jd [u]},  then  d[x\  =  d'[x\  and 
d'[x]  =  min„ee{d'[t;]}. 

Proof:  Proof  by  contradiction.  Assume  that  d'[x\  <  d[x], 
then  x  has  to  use  e  on  its  new  shortest  hyperpath.  Since  we 
consider  only  simple  hyperpaths  and  x  £  e,  we  have  E[x]  =  e. 
Therefore  its  parent  y  =  P\v]  cannot  use  e  on  its  shortest 
hyperpath,  which  implies  that  the  shortest  distance  to  y  does 
not  change:  d[y]  =  d'[y].  Given  that  y  is  the  parent  of  x  on  its 
new  shortest  hyperpath,  we  have  d[y]  =  d'[y]  <  d'[ x]  <  c?[a;] 
which  contradicts  to  the  definition  of  x. 

For  the  second  statement,  assume  there  exists  z  £  e  such 
that  d'[z)  <  d'[x\.  Based  on  the  definition  of  x  and  the 
hypothetical  assumption,  d[z\  >  d[x\  =  d'[x\  >  d'\z).  It  thus 
follows  that  z’s  shortest  hyperpath  changes  and  E\z\  =  e  in  the 
new  shortest  hyperpath.  Follow  the  same  line  of  arguments  by 
considering  the  parent  of  z,  we  arrive  at  the  same  contradiction 
in  terms  of  the  definition  of  x.  ■ 

Lemma  3:  For  any  vertex  v,  v  is  enqueued  into  Q  if  and 
only  if  d'[v]  <  d[v]. 

Proof:  Consider  first  that  v  is  enqueued  into  Q.  From  the 
algorithm,  this  can  only  happen  if  there  exists  a  neighbor  z 
and  a  hyperedge  e  9  v,  z  such  that  D[z\  +  w(e)  <  D[v\.  We 
thus  have  d[i>]  >  D[v\  >  D[z\  +  w(e)  >  d'[v\  (note  that  at 
any  time,  d[u]  >  D[v\  >  d'[v\,  which  can  be  easily  seen  from 
the  procedure  of  the  algorithm). 

We  now  prove  the  converse.  Assume  that  d'[v\  <  d[v}. 
Let  p  =  {ei,  e2, . . . ,  ei,  e,  e*+i, . . . ,  e;}  be  v’s  new  shortest 
hyperpath.  There  exists  14+1  &  e  FI  ei+i  such  that  d'[iii- )_i]  < 
d[ui+ 1].  In  Step  1  of  the  algorithm,  u.l+\  is  enqueued.  Simi¬ 
larly,  there  exists  m+2  £  ej+i  Clei+2  with  d'[ui+ 2]  <  d[ui+ 2]. 
Then  itj+ 2  will  be  enqueued  in  Step  2  of  the  algorithm  when 
Ui+ 1  is  dequeued  if  it  has  not  been  enqueued  before  that. 
Repeating  this  line  of  argument,  we  conclude  that  there  exits 
ui  £  ei- 1  FI  ei  with  d'[ui]  <  d[ui]  and  ui  is  enqueued  into 
Q.  Then  v  will  be  enqueued  when  iq  is  dequeued  if  it  is  not 
enqueued  already.  ■ 

Lemma  4:  For  each  v  dequeued  from  Q,  D[v ]  =  d'[v\. 

Proof:  We  first  show  that  if  u  is  dequeued  before  v,  then 
D[u ]  <  D[v]  at  the  instants  when  they  are  dequeued.  We  prove 
this  by  induction.  The  initial  condition  holds  trivially.  Then 
assume  it  is  true  for  the  first  l  dequeued  vertices  ...  ,zi. 
Consider  the  ( l  +  l)th  dequeued  vertex  zi+ 1.  At  the  instant 
when  zi  is  dequeued,  if  D[zi+ 1]  is  updated  based  on  D[zi ] 
in  Step  2,  then  D[zi ]  <  D[zi+ 1]  even  after  the  update.  If, 
on  the  other  hand,  D[zi+ 1]  is  not  updated  at  this  instant, 
then  D[zi]  <  D[zi+ 1]  given  that  the  dequeued  vertex  has  the 
smallest  distance. 

Next,  we  prove  the  lemma  by  induction.  From  Step  1  of  the 
algorithm,  all  the  affected  vertices  v  in  e  will  be  dequeued  first 
with  E[v]  =  e,  P[u]  =  x,  and  D[v]  =  d'[x\  +  w{e).  Based  on 
Lemma  2,  D[v\  <  d'[u\  +  w(e)  for  any  u  £  e.  It  thus  follows 
that  the  hyperpath  to  v  through  x  and  e  is  the  shortest  one 
with  D[v]  =  d'[vj. 

Assume  for  z\,...,zi,  D[zi\  =  d'[zi\  are  satisfied  for  all 
i  =  1, . . .  ,1.  Consider  the  (/  +  l)th  dequeued  vertex  zi+ 1  ^  e. 
Let  u  =  P[zi+ 1]  be  its  parent  in  the  new  shortest  hyperpath. 


Then  based  on  the  fact  that  distances  of  the  dequeued  vertices 
are  monotonically  increasing  with  the  order  of  the  dequeueing 
as  shown  at  the  beginning  of  the  proof,  u  cannot  be  any 
vertex  dequeued  after  zi+  \ .  Since  zi+ 1  (f  e,  it  is  also  clear 
that  u  cannot  be  an  unaffected  vertex(otherwise,  z/+]  will 
be  unaffected,  which  contradicts  Lemma  3).  We  thus  have 
u  £  {z±, . . . ,  zi}.  Let  u  =  Zi.  Then  when  Zi  is  dequeued  , 
D[zi- t_i]  will  be  updated  to  the  shortest  distance  d'\zi+ 1]  due 
to  the  induction  hypothesis  of  D[zi]  =  d'\zi}.  This  completes 
the  proof.  ■ 

Based  on  Lemma  3  and  4,  the  shortest  distances  of  all  affected 
vertices  will  be  updated  correctly.  Based  on  Lemma  3,  all 
unaffected  vertices  will  not  be  enqueued,  and  their  distances 
remain  the  same.  It  is  not  difficulty  to  see  from  the  algorithm 
that  P[u]  and  E\v]  are  also  correctly  maintained  for  all  v. 

Appendix  B:  Proof  of  Theorem  2 

We  first  show  the  correctness  of  the  coloring  process  as 
given  in  the  following  lemma. 

Lemma  5:  The  coloring  process  correctly  colors  all  the 
affected  vertices. 

Proof:  We  first  state  the  following  simple  facts  without 
proof:  given  a  relationship  tree,  after  the  hyperedge  weight 
increase,  (1)  if  v  is  pink  or  white,  then  all  its  descendent  in 
this  relationship  tree  are  white;  (2)  if  a  v  is  red,  then  all  its 
children  in  the  relationship  tree  are  either  pink  or  red;  (3)  if  a 
v  is  affected,  either  v  £  e  or  P[v\  is  red.  These  facts  can  be 
directly  obtained  from  the  definition  of  the  color.  It  is  also  easy 
to  see  that  vertices  are  dequeued  from  M  in  a  nondecreasing 
order  of  their  current  distance  D[-}.  This  is  because  each  time 
a  vertex  z  is  dequeued  from  M,  the  possible  new  vertices  to 
be  enqueued  into  M  are  z’s  children  with  distances  no  smaller 
than  D[z}. 

Then,  the  proof  of  the  lemma  has  two  parts:  first  we  prove 
that  all  affected  vertices  are  enqueued  into  M;  then  we  prove 
by  induction  that  only  affected  vertices  are  enqueued  into  M 
and  their  colors  are  correctly  identified. 

We  prove  the  first  part  by  contradiction.  Assume  that  there 
exists  an  affected  vertex  v  that  is  not  enqueued  into  M.  It  is 
easy  to  see  that  v  ^  e  because  all  the  affected  vertices  in  e 
are  enqueued  in  Step  1.  Based  on  the  third  fact  stated  above, 
P[v]  is  red.  Based  on  the  hypothesis,  P[u]  is  not  enqueued 
(otherwise,  v  will  be  enqueued  in  Step  2).  Continue  this  line 
of  arguments,  we  eventually  reach  the  root  of  the  relationship 
tree  and  arrive  at  the  contradiction  that  the  source  s  is  red. 

We  prove  the  second  part  by  induction.  It  is  easy  to  see  that 
all  the  vertices  initially  enqueued  into  M  are  affected  vertices. 
It  remains  to  show  that  the  first  vertex  Z\  dequeued  from 
M  is  colored  (pink  or  red)  correctly.  To  show  that,  we  need 
to  establish  that  the  algorithm  correctly  determines  whether 
there  is  an  alternative  shortest  hyperpath  to  z\  with  the  same 
distance,  i.e.,  d[z{\  =  d'[zi].  The  key  here  is  to  show  that 
checking  the  currently  non-red  neighbors  (which  may  become 
red  in  the  future)  of  z\  will  not  lead  to  a  false  alternative  path. 
This  follows  from  the  fact  that  z\  has  the  smallest  distance  D[-] 
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among  all  affected  vertices  (which  belong  to  the  set  of  vertices 
consisting  of  the  affected  vertices  in  e  and  their  descendents). 

Next,  assume  that  vertices  zi,  Z2, . .  ■ ,  zi  dequeued  from  M 
are  all  affected  vertices  and  are  correctly  colored.  Consider  the 
next  dequeued  vertex  zi+ 1.  It  is  an  affected  vertex  because  it  is 
either  enqueued  in  Step  1  with  E[v\  =  e  or  enqueued  in  Step  2 
with  a  red  parent.  To  show  that  z/+  \  will  be  colored  correctly, 
we  use  a  similar  argument  by  showing  that  the  currently  non- 
red  neighbors  of  zi+ 1  will  not  give  a  false  alternative  path. 
The  latter  follows  from  the  fact  that  all  affected  vertices  will 
be  enqueued  and  those  dequeued  after  zi+\  have  distances  no 
smaller  than  D[zi+ 1].  This  completes  the  induction.  ■ 

We  now  show  that  D[v\,  P[v]  and  E[v]  are  correctly 
maintained  for  all  v.  For  each  red  vertex  v,  its  distance  is 
set  based  on  the  current  shortest  distance  from  a  non-red 
neighbor  in  Step  3. a.  The  rest  of  the  algorithm  is  essentially 
Gallo’s  extension  of  Dijkastra’s  algorithm  with  the  current 
initial  distance.  The  correctness  of  the  algorithm  thus  follows. 
It  is  not  difficult  to  see  that  P[-]  and  E[-\  are  correctly  updated 
for  both  red  and  pink  vertices. 
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